Autonomous Stochastic Resonance in Fully Frustrated Josephson-Junction Ladders 
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We investigate autonomous stochastic resonance in fully frustrated Josephson-junction ladders, 
which are driven by uniform constant currents. At zero temperature large currents induce oscillations 
between the two ground states, while for small currents the lattice potential forces the system to 
remain in one of the two states. At finite temperatures, on the other hand, oscillations between the 
two states develop even below the critical current; the signal-to-noise ratio is found to display array- 
enhanced stochastic resonance. It is suggested that such behavior may be observed experimentally 
through the measurement of the staggered voltage. 
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The stochastic resonance (SR), which describes the 
phenomenon that a noise of finite strength optimizes the 
response signal to a weak periodic external force in a 
nonhnear systemj^has attracted much attention for the 
past two decades. □ Conventional SR phenomena are gen- 
erally known to require an energy barrier between two 
stable states, a weak external periodic force, and noise. 
However, the SR behavior was also exhibited wiijhout 
any external periodic force, in limit-cycle systemsH and 
excitable systems.cru Subsequently such coherence reso- 
nance or autonomous SR has also been observed in sys- 
tems with dela}cl and inertia.tl In the latter cases, SR 
phenomena can still be understood by the time-scale 
matching argument since the delay or the inertia pro- 
vides an external time-scale. In the former, on the other 
hand, such an external time scale apparently does not 
exist, and different noise dependencies of the activation 
time and the excursion time have been proposed as an 
explanation of the coherence resonance.Q Recently, the 
fully-frustrated Josephson-junction ladder has been pro- 
posed to give a good physical realization of the standard 
two-state SR system with many degrees of freedom.Q In 
the presence of the external current which is periodic in 
time and staggered in space, the SR behavior and other 
rich physics have been demonstrated theoretically in the 
ladder system. Furthermore, recent progress in the fabri- 
cation technique makes such Josenhson-junction systems 
available for experimental study.ErLi 

In this work we study fully frustrated Josephson- 
junction ladders driven by uniform constant currents, 
paying particular attention to the possibility of the au- 
tonomous SR phenomena. Note that this system is 
bistable, possessing two stable states; this is in sharp con- 
trast with the usual excitable system, characterized by a 
single stable state. Numerical integration of the coupled 
equations of motion, established within the resistively 
shunted junction (RS J) model, shows the existence of the 
critical current beyond which oscillations emerge between 
the two ground states. Below the critical current, such 
oscillations are suppressed by the lattice potential; here 



the addition of noise currents, relevant at finite temper- 
atures, induces again oscillations, giving rise to SR be- 
havior. In particular, the signal-to-noise ratio (SNR) is 
found to display array-enhanced SR phenomena, which 
may be observed experimentally through the measure- 
ment of the staggered voltage. To our knowledge, such 
array-enhanced SR behavior has not been observed in 
bistable system without periodic driving although coher- 
ence resonance was reported in the bistable regime of the 
FitzHugh-Nagumo model.EJ 

The RSJ model, with single-junction critical current Ic 
and shunt resistance R, is described by the set of equa- 
tions of motion for the phaSj(ji,0i of the superconducting 
order parameter on grain ii 
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where (f>ij = (pi — (f)j is the phase difference across the 
junction (ij), the thermal noise current rjij is assumed to 
be the white noise satisfying 
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S{T){S,kSji ~ SaSjk), (2) 



and the primed summation runs over the nearest neigh- 
bors of grain i. We choose the Landau gauge and write 
the gauge field Aij, which describes the transverse mag- 
netic field, in the form 



Ai-j — 



for {ij) legs, 
27: fx for (ij) rungs. 
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where / = ^/^o is the flux <& per plaquette in units of the 
flux quantum $o = hc/2e and x {= 1,2,... , L) denotes 
the position along the leg. As illustrated in Fig. the 
ladder, placed along the x direction, is driven by the 
uniform constant currents / along the rung direction (i.e., 
the y direction): The external current /™* fed into grain 
I is thus given by 
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FIG. 1: Schematic diagram of a fully frustrated Josephson 
junction ladder driven by uniform constant currents. 
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FIG. 2: Time evolution of the staggered magnetization for 
various currents / = 0.2, 0.25, and 0.3 from above. For clarity, 
the data corresponding to 7 = 0.2 and 0.3 have been shifted 
by 1.5 upward and downward, respectively. 



where y {= 1,2) is the leg index of the ith grain. Hence- 
forth, we write the current, the temperature, and the 
time in units of Ic, hlcf^e, and h/2eRIc, respectively. 

In the absence of external currents, the fully frustrated 
(/ = 1/2) ladder is well known to have two degenerate 
ground states, where a vortex is located on every two 
plaquettesE (Note that the length L of a fully frustrated 
ladder should be an even number.) These states can be 
characterized by the staggered magnetization 
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where nx is 
vorticity is 
(1/2^) Ep, 



the vorticity 
given by the 



on plaquette X. The 
plaquette sum nx = 
1/2, where the gauge- 



invariant phase difference is defined modulo 2tt in the 
range (— tt, tt]. Here uniformly applied currents in the 
y direction drive vortices to the x direction, generating 
a flow of the vortex array along the ladder. Due to the 
spatial periodicity of the vortex array, the flow in the x di- 
rection results in the alternation between the two ground 



states, which can be manifested by the oscillation of the 
staggered magnetization. In contrast to a homogeneous 
thin film, the array system has a lattice potential which 
tends to suppress the flow of a vortex array. Accordingly, 
oscillations of the staggered magnetization should appear 
only when the external currents exceed the critical value 
Iq = y^—2.t3 The time evolution of the staggered mag- 
netization at zero temperature, obtained from direct in- 
tegration of Eq. (|l|) and presented in Fig. || for various 
values of /, indeed confirms the existence of the critical 
value: For small currents (/ = 0.2 < Jq), the staggered 
magnetization remains constant, indicating that the lad- 
der is pinned by the lattice potential to one of the two 
ground states. Above the critical value /o, there arise os- 
cillations of the staggered magnetization. Note also that 
the oscillation frequency increases with currents, reflect- 
ing that the vortex array flows faster at larger currents. 

To probe the possibility of autonomous SR here, we 
consider the average rate of the change in the staggered 
magnetization 
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which measures the transition rate between the ground 
states. The presence of periodic oscillations between the 
ground states is manifested by the peak in the power 
spectrum Sy{uj) of v{t) at the corresponding frequency. 
We then compute the SNR: 



i? EE 10 logi 
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where the signal S = Sy{uj—ujp) is the peak value of the 
power spectrum (at the peak frequency LOp). The back- 
ground noise level N is estimated by the power spec- 
trum averaged over the interval of frequency much higher 
than Wp, where the power spectrum is flat. Equation (|l|) 
has thus been integrated directly with the time steps 
Ai = 0.05, from which v{t) has been computed. We 
have also reduced At to 0.01, to obtain the same results 
within the numerical accuracy, and followed an annealing 
schedule with the equilibration time 1000 at each tem- 
perature. The data have been averaged over 200 to 5000 
independent runs, depending on the system size. For 
convenience, 5t has been set to be 0.25. We have also 
considered smaller values of (5t, and found that the over- 
all results remain qualitatively the same although higher 
noise levels make it more formidable to perform precise 
numerical analysis. The periodic boundary conditions 
have been imposed along the x direction. 

We first examine the effects of noise above the critical 
current. Figure ^ shows the power spectrum at various 
temperatures in the ladder of length L = 50 above the 
critical current (/ — 0.25). The power spectrum exhibits 
a sharp peak even at zero temperature, arising from the 
current-induced oscillation. As the temperature is raised, 
however, the signal S decreases monotonically, reflecting 
that the thermal noise disturbs the coherent motion of 
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FIG. 3: Power spectrum of the transition rate v{t) (in 
arbitrary units) in a ladder of length L = 50 for I — 
0.25 at various temperatures T = 0.01, 0.02, 0.03, 0.04, 
0.06, 0.08, 0.12, 0.16, 0.24, and 0.4 from left. 
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FIG. 4: SNR versus temperature in a ladder of length L = 50 
for / — 0.25. Here and in other figures applicable, error bars 
have been estimated by the standard deviation, and the lines 
connecting data points are merely guides to eyes. 

the vortex array. Accordingly, the SNR displays mono- 
tonic decrease with the temperature, as shown in Fig. ^. 

We next reduce the current / to 0.2, below /q. Unlike 
at zero temperature, there emerge oscillations between 
the two ground states at finite temperatures, as man- 
ifested by the peak of Sy{uj) in Fig. ^. The signal S 
increases as the temperature is raised from zero while 
further increase of the temperature tends to suppress the 
peak, eventually forcing the power spectrum to be white. 
The SNR plotted in Fig. | for / = 0.2 and lower values 
clearly displajwSi-iautonomous SR, i.e., SR without peri- 
odic driving.ffl'Q Here it should be stressed that the SR 
observed above does not describe the direct response to 
the external dc current. As shown in Fig. ^, the peak fre- 
quency, which is proportional to the mean flow velocity of 
the vortex array, increases monotonically with the tem- 
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FIG. 5: Power spectrum of v(t) (in arbitrary units) in a lad- 
der of length L = 50 for I — 0.2 at various temperatures 
T = 0.01, 0.02, 0.03, 0.04, 0.06, 0.08, 0.12, 0.16, 0.28, and 0.4 
from below. 
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FIG. 6: SNR versus temperature in a ladder of length L = 50 
for various currents below Iq. 

perature, exhibiting no SR behavior. This is consistent 
with the previous study which concluded that the mobil- 
ity of an overdamped particle-4p a washboard potential 
does not display SR behavior 113 

Similar SR phenomena are also observed in the behav- 
ior of the staggered magnetization m(t) itself. Figure ^ 
displays the power spectrum Sm{uj) of the staggered mag- 
netization m{t) both above and below the critical cur- 
rent. Above the critical current, indeed similarly to Sy (uj) 
shown in Fig. ^ S',„(ci;) exhibits a sharp zero-temperature 
peak at a finite frequency, which gets broadened gradu- 
ally with the introduction of thermal noise. Below the 
critical current, on the other hand, the zero-temperature 
peak develops at zero frequency {uj — 0), indicating the 
absence of oscillation in the staggered magnetization. As 
the temperature is raised from zero, the dc component 
Sm{'-o=0) becomes smaller; instead there emerges a peak 
at a finite frequency ujp, which, like the peak frequency 
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FIG. 7; Peak frequency ujp as a function of the temperature 
in a ladder of length L = 50 for / = 0.2. 



ujp of Sv{u}) (see Fig. R), depends on the temperature. 
These behaviors of Smluj) and those of are in gen- 

eral consistent with the relation Si,{Ld) cx uo'^Smi^), accu- 
rate for small 5t. Note here that such development of a 
finite- frequency peak at finite temperatures is an indica- 
tion of the coherence resonance phenomena observed, in 
the bistable regime of the Fitzhugh-Nagumo model.EHl 

To quantify conveniently the expected SR, behavior of 
m{t), we compute the coherence measureB'Ej 



/3 = Smi^^p) ( ^ 



(8) 



-0J„ 



where the half- width of the peak is given by Aw = 
with LOr (> '^p) satisfying 25'm(<^r) = Sra{^p)- The co- 
herence measure (3, which is plotted as a function of 
temperature in Fig. ^, also reveals the nature of the au- 
tonomous SR behavior in the system. Above the critical 
current it decreases exponentially with the temperature 
[see Fig. ^(a) for / — 0.25] while below the critical cur- 
rent a prominent peak appears at a finite temperature, 
as demonstrated in Fig. ||(b) for / = 0.2. Thus confirmed 
is the autonomous SR behavior of the staggered magne- 
tization m{t). 

The autonomous SR in the fully frustrated ladder is 
rather peculiar, arising from the combination of the ex- 
ternal uniform currents and the spatial periodicity of the 
ground state inherent in the system. Although currents 
smaller than /q do not induce the vortex array to move, 
they add an overall gradient to the lattice potential along 
the X direction. At finite temperatures, the additional 
noise currents assist vortices to hop to neighboring pla- 
quettes while the potential gradient generated by the ex- 
ternal currents makes hopping in one direction more fa- 
vorable than that in the other. Consequently, the vortex 
array is encouraged to flow to one direction, leading to 
oscillations of the transition rate v{t). These oscillations 
grow with the temperature, which gives the enhancement 
of the signal. At higher temperatures, however, the lat- 
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FIG. 8: Power spectrum of the staggered magneti- 
zation m{t) (in arbitrary units) at various temperatures 
T = 0.01, 0.02, 0.03, 0.04,' 0.06, 0.08, 0.12, 0.16, 0.24, and 0.4 
(from left) in a ladder of length L — 50 for (a) / — 0.25 and 
(b) / = 0.2. 



tice potential gradient can be neglected, and the random 
nature of the noise disturbs the coherent flow, weaken- 
ing the signal. In this manner the two different roles pf 
the noise produces the autonomous SR in the system,0 
which also provides a natural explanation of the current- 
dependence of the SNR behavior in Fig. ^: At small cur- 
rents the lattice potential barrier between neighboring 
sites remain rather high, and the noise required to over- 
come the barrier should be strong. Accordingly, as the 
applied currents are reduced, the SR temperature is ex- 
pected to increase while the SNR peak diminishes due to 
the large incoherent contributions of the noise. 

Figure showing the SNR versus temperature for 
various lengths, discloses another interesting feature of 
the SR behavior in the ladder system. In the ladder of 
length L — 2, which is the minimum length required 
for a fully frustrated ladder, the SNR exhibits no sign 
of SR. As the length is increased, there first appears a 
broad peak in the SNR, which then develops into a promi- 
nent peak. Such array- enhanced SR phenomena reflect 
the stronger coherence between vortices in longer lad- 
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FIG. 9: Coherence measure versus temperature in a ladder of 
length L = 50 for (a) / = 0.25 and (b) I = 0.2. 
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FIG. 10; SNR versus temperature for various sizes and / 
0.2. 




FIG. 11: SNR for the staggered voltage as a function of the 
temperature in a ladder of length L = 50 for I = 0.2. 



ders. However, the enhancement should not persist with 
the length since the ladder does not evolve long-range 
order in the thermodynamic limit. Indeed the SNR in 
Fig. ^ eventually saturates for L > 30, appearing inde- 
pendent of the length. While array-enhanced resonance 
was first renorted in the array systems exhibiting conven- 
tional SRJIj recent studies have also revealed the possi- 
bility of array-enhanced coherence resonance in coupled 
excitable systems.EJ It is remarkable that similar array- 
enhancement of autonomous SR can also be observed in 
the Josephson-junction ladder, which possesses two sta- 
ble states like the conventional SR system (but in the 
absence of periodic driving). 

Finally, we discuss how to observe the autonomous SR 
in experiment. As a candidate for the measurable quan- 
tity characterizing the autonomous SR, we suggest the 
average staggered voltage: 



Vs{t) 



(9) 



where Vx{t) is the rung voltage at position x averaged 
over the interval 5t around time t. The rigid motion of 
the vortex array at low temperatures should give rise to 
periodic oscillations of Vs{t), similarly to the transition 
rate v{t). The SNR for the staggered voltage Vs{t), de- 
noted by Rs, is plotted in Fig. where SR behavior 
similar to that of v{t) is shown. We thus suggest that 
the autonomous SR phenomena in the ladder should be 
observed through the measurement of the staggered volt- 
age. 

In summary we have investigated fully frustrated 
Josephson-junction ladders, driven by uniform constant 
currents, with regard to the possibility of stochastic res- 
onance. While the system under large currents displays 
oscillations between the two ground states, the lattice 
potential suppresses such oscillations for small currents. 
Still the addition of noise currents, relevant at finite 
temperatures, induces again oscillations, giving rise to 
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stochastic resonance behavior. In particular the signal- 
to-noise ratio has been found to display array-enhanced 
stochastic resonance phenomena. It has also been sug- 
gested that such behavior may be observed experimen- 
tally through the measurement of the staggered voltage. 
In view of this, the fully frustrated Josephson-junction 
ladder makes a good physical realization of the system 
with many degrees of freedom, adequate for the study 



of the autonomous as well as the conventional stochastic 
resonance phenomena. 
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